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of  the  lossy  half-space  on  their  input  impedance  and  the  far-field  radiation 
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1.  INTRODUCTION 

The  conventional  approach  to  analyzing  antenna  structures  radiating 
in  the  presence  of  a lossy  half-space  involves  repeated  evaluation  of  the 
Sommerfeld  integrals,  which  were  originally  introduced  by  Sommerfeld  about 
70  years  ago  [1]  and  appeared  in  the  expressions  for  the  vector  potentials  [2]. 
Since  these  infinite  integrals  are  generally  highly  oscillatory  and  difficult 
co  evaluate  numerically,  much  attention  has  been  focused  in  recent  years 
on  developing  techniques  for  efficiently  evaluating  the  Sommerfeld 
integrals  without  unduly  sacrificing  the  accuracy  [3]  - [14].  Even  though 
the  latest  reported  procedures  [7]  - [12]  require  an  order  of  magnitude  less 
computing  time  than  the  earlier  "brute-force"  numerical  integration  tech- 
niques [3]  - [5],  the  overall  computing  time  severely  limits  the  physical 
dimensions  of  the  antenna  structures  being  numerically  analyzed. 

Brittingham  et  al.  [10]  have  used  an  interpolation  scheme  on  a precalculated 
grid  of  Sommerfeld  integral  values  in  order  to  analyze  larger  structures. 
However,  this  technique  is  only  useful  for  fixed  frequency  and  ground 
parameters,  since  a new  grid  is  required  each  time  any  of  these  parameters 
are  changed. 

Considerable  computing  time  is  saved  if  one  uses  the  first  term  in  the 
asymptotic  expansion  of  the  Sommerfeld  integrals,  better  known  as  the 
Fresnel's  Reflection  Coefficient  Method  (RCM) . These  approximations, 
which  are  expressed  in  a simple  closed  form,  are  valid  only  when  the  antenna 


tfl 

> 


i 


J 


structures  and  the  observation  points  are  sufficiently  high  above  the  ground. 
As  expected,  the  RCM  expressions  have  only  been  employed  for  analyzing 
antenna  structures  at  the  high  end  of  the  frequency  spectrum  [10],  [12],  [13]. 


An  important  contribution  of  this  paper  is  tbe  development  of  a novel 
technique  which  is  computationally  comparable  to  the  ROM  approximation, 
yet  valid  for  a much  wider  range  of  parameters. 

In  Chapter  2,  the  complex  vector  potential  expressions  for  an 
arbitrarily  oriented  electric-current  element  source  over  a lossy 
half-space  are  derived  demonstrating  that  three  of  the  vector  potential 
components  contain  the  troublesome  Sommerfeld  Integrals.  An  efficient 
numerical  technique,  based  on  the  Steepest  Descent  Path  (SDP)  integration, 
is  presented  in  Chapter  3 for  evaluating  the  exact  Sommerfeld  integral 
expressions,  and  the  results  are  compared  with  the  one-  and  two-term 
asymptotic  expansions  of  these  integrals.  In  Chapter  4,  a unique 
procedure  is  developed  in  which  the  well-behaved  Fourier  transform 
representations  of  the  Sommerfeld  integrals  are  approximated  such  that 
the  inverse  transform  is  performed  via  a set  of  known  exact  identities. 

These  approximate  expressions  are  shown  not  only  to  closely  follow  the 
SDP  integration  results  for  a wide  range  of  parameters  of  practical  interest, 
but  to  be  in  a convenient  form  for  numerical  evaluation.  Finally,  in 
Chapter  5,  several  examples,  which  are  solved  via  the  method  of  moments 
in  conjunction  with  the  approximate  expressions  derived  in  Chapter  -t,  are 
presented  of  various  antenna  structures  radiating  over  a lossv  half-space. 
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ARBITRARY  CURRENT  ELEMENT  RADIATING  OVER  LOSSY  HALF-SPACE 


The  geometry  of  ;in  arbitrarily  oriented  current  element  P ^ over  .1 


lossy  half-space  is  depicted  in  Figure  1.  Regions  l and  2 are  characterized 
by  (Cj  - * Uq)  and  (e->  “ e:>re()’  u->  " 1 V ’ respectively,  where  tQ 


and  u(1  are  free-space  parameters.  In  addition  to  the  standard  Cartesian 


coordinate  svstem  (x,  v,  z),  two  spherical  coordinate  systems  (r,  0^,  |^) 


and  (r,,  4>.,)  are  also  defined  in  Figure  1 centered  about  the  source 


point  Pj  and  its  geometrical  image  point  P,,  respectively.  Our  objective 


is  determine  the  field  radiated  by  P^  at  the  observation  point  0,  in 


the  presence  of  the  lossy  half-space  (region  2). 

2.  I Vector  Potential  Approach 

Starting  with  Maxwell's  equation  and  tin’  suppressed  time  convention 
exp(jiot),  viz.. 


V x H » jioe_c  E + .1 
Or 


(2.1a) 


V x E » -jiou^H 


(2.  lb) 


one  may  define  the  vector  potential  H as 


H - J<0t'  £ V x FI 
J Or 


(2.2) 


Introduction  of  a scalar  potential  i'  from 


V4>  « E - 10^11, 11 
0 1)  r 


(2.3) 


and  application  of  the  Lorentz  gauge 

v • n - - 0 

allows  one  to  finally  express  Maxwell's  equation  as 


(2.s) 


(V  + ion  - -d<ot- c ) * .1 

Or 


(2.  ') 
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H = jtoe,,t:  V x n (2.6) 

J Or 

E - (W*  + k2)IT  (2.7) 

7 •) 

where  k~  » . The  preceding  results  are  general  and  valid  for  both 

regions  l and  2.  The  boundary  conditions  needed  to  solve  the  vector  dif- 
ferential Equation  (2.5)  can  be  obtained  simply  by  enforcing  the  continuity 
of  the  tangential  E-  and  H-field  components  at  the  interface. 

As  in  most  inf inite- inter face-type  problems  reported  in  the  literature, 
the  Fourier  transform  technique  is  employed  for  deriving  the  vector 
potential  expressions.  The  two-dimensional  Fourier  transform  pair  is 
defined  as 


if 


JO 

J"  a exp [-j  (ax  + By)  ] dx  dy 

— OU 


(2.8a) 


-•IT 


i_  r I 

- i J 


exp[J(ax  + By) 1 da  dS 


(2.8b) 


Throughout  this  papet,  ' on  top  denotes  the  transformed  quantity. 

Without  loss  of  general  itv,  the  problem  of  an  arbitrarily  oriented 
current  source  over  lossy  half-space  (Figure  1)  can  be  treated  as  two 
independent  cases:  first,  the  problem  of  a vertical  current  source 
oriented  in  the  2-direction  and, second,  a horizontal  source  in  the 
x-directlon.  These  two  cases  are  analyzed  in  detail  in  the  following 
sections  with  subscripts  v and  h denoting  the  field  quantities  belonging 
to  the  vertical  and  the  horizontal  current  sources,  respect ive lv . Once 
the  vector  potential  expressions  are  derived,  the  E-  and  H-field  components 
are  directly  obtained  bv  transforming  the  general  vector  equations  given 
in  (2.6)  and  (2.7)  into  the  following  convenient  matrix  forms: 


, , 


...  ....  


( 


1 

k“  - a“ 


-at3  k“  - S' 


9 

Ja  3l 


Ja9l  nx 


jBk  nv 


<k2 + y l3> 


0 

1 N 

l 

J3 

3_ 

3z 

0 

-ja 

-JB 

Jot 

0 

(2.10) 


2 . 2 Vertical  Current  Element -Transform  Domain  Expressions 

The  current  element  Is  assumed  to  be  in  the  z-direction  (see  Figure  1) 
located  at  height  h over  the  lossy  ground,  and  with  a moment  I^dz': 


J - z I dz'  S(x)  <5(v)  <S(z) 

V V 


(2.11) 


As  has  been  verified  in  the  literature  [ 2 ] > only  the  vertical  component 
of  the  vector  potential  is  needed  to  satisfy  all  of  the  boundary 
conditions,  namelv 


n - ill  ; i - 1,2 

vi  viz 


(2.12) 


where  i « 1,2  indicates  the  region  under  consideration.  Using  the  source 

condition  in  (2.11),  one  can  express  the  general  solution  for  the  vector 

potential  FI  . bv  satisfying  both  the  wave  Equation  (2.5)  and  the 
viz 


radiation  conditions  as 

^vlz  3 Iv0  exp(-jY12) 


0 (2.13a) 


H , * A,  exp( jv ,z)  ; 2^0 


(,  2 . I 3b) 
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where  A^  and  A.,  are  arbitrary  constants  and 
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(2.14) 

(2. 15) 
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i - 1.2  ( 2 . In) 

It  the  continuity  of  the  tangential  components  of  the  E-  and  H-tields  across 


the  boundary  is  enforced,  the  following  constraints  on  FI  are  obtained 


f n . I - t-  n , | 

a 2 v l * 3 z v 2 z 1 

1 'z-0  'z-0 
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viz1 
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(2. 1 7h) 
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where  the  complex  constant  < is  defined  as 


< * e /c 

2r  Ir 


(2. 18) 


Constants  A^  and  A.,,  present  in  the  vector  potential  expressions  (2.13a)  and 


(2.13b)  can  be  determined  by  satisfying  the  boundary  conditions  (2.17a)  and 
(2.17b) : 
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t - )T 
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vO  J(k.A,+V,) 
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l-  -J 

exp  (-h  jh) 
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If  (2. 1'))  is  substituted  back  into  (2.13a)  and  (2.1!b),  the  complete  vector 
potential  expression  in  the  transform  domain  can  be  written  in  the  following 
form: 


i r 

n - it  + t + n 

viz  viz  v 1 z 0 v 1 : 


(2.20) 


where 


„i 


(2.21a) 


I exp  | - j n , \z  - hi  ] /2 


nvU  “ " lvO  expl'jYl(Z  + h)  I / 2 j Y j 

onvi2  ’ lyO  jlTy'+vT)  expl~jYl(z  + h) 


(2. Jib) 
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nv22  ■ [V()  + y"‘)  exp  (-jYlh)  exp  (i>22)  • (2-22) 

Note  that  the  first  two  terms  of  II  , , viz.,  II*  and  H1”,  , can  be 

viz  viz  viz 

interpreted  as  the  direct  and  the  perfect  ground  reflection  contributions 
to  the  fields  for  observation  points  in  region  1.  Thus,  the  remaining 
component  ^JI  ^ Is  simply  the  correction  term  to  the  perfect  ground 
solution  in  the  general  lossy  half-space  problem. 

2 . 3 Horizontal  Current  Element-Transform  Domain  Expressions 

The  current  element  is  assumed  to  be  in  the  x-direction  (see  Figure  1), 


located  at  height  h over  lossy  ground  and  with  a moment  I^dx’ 
* x I^dx’  »S(x)  d'(y)  3(z  - h) 


l 2. 2 3) 


As  expected,  the  horizontal  current  element  solution  above  a lossy  half- 
space is  indeed  more  complicated  than  that  for  the  vertical  case,  and 
two  vector  potential  components  are  needed  to  obtain  a complete 


solution  [2],  namely, 

i - x n. . + z n. . ; i - 1,2 

hi  hix  hiz 

Using  the  source  condition  in  (2.23),  the  general  solutions  for  the 


(2.2s) 


vector  potential  components  satisfying  both  the  wave  Equation  (2.5)  and 
the  radiation  condition  can  be  expressed  as 


1,,0  *xp[-j  , j ]z  - h | 1 / 2 J ) 


•xpi-h  .=)  : •'  1 11 
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exp  ( j Y , z ) ; z < 0 


where  >,  is  Jet iued  in  (2.15)  and  1,  . is  assumed  to  he 
i hO 


(2.  26) 


I * ijtot'  € ) I dx' 
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If  the  continuity  of  the  tangential  components  of  the  E-  and  the  H-t'ields 

across  the  boundary  is  enforced,  the  following  constraints  on  H,  , and 

hix 

FI,  , are  obtained: 
hiz 


JOt  3^  "Ja 


1 0 

0 1 


o n. 
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The  conditions  in  (.2.2$)  are  used  to  determine  the  constants  3.  and  3 

tx  ic 

present  in  the  general  vector  potential  expressions  given  in  i2.25) 
and  (2.26): 
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If  (2.29a)  and  (2.29b)  are  substituted  back  into  (2.25)  and  (2.26),  the 
complete  transform  domain  vector  potential  expressions  can  be  put  into 
the  following  convenient  form: 


n, . = nf,  + rt,r.  + -IT. . 

hlx  hlx  hlx  0 hlx 


(2.30) 


where 


Tilx  = rh0  exp  [“jYll2  - hll/2jY1 


(2.31a) 


nhlx  = ' Th0  exp  ["jYl(z  + h)]/2jYl 


(2.31b) 


0nhlx  = Ih0  !(—+—)  eXP  ['jYl(Z  + h)1 


(2.31c) 
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(2.32) 


nh2x  = Xh0  jV(77+72)  eXP  (“jYlh)  eXP  (jY2Z) 

Y1  ' Y2 

\2z  = rh0  ja  72~ exp  (-jYlh)  exp  (jY2z) 

k1K(KY1  + y2) 


(2.33) 


(2.34) 


As  in  the  vertical  case,  the  first  two  terms  of  the  dominant  vector 

i it 

potential  component,  viz.,  and  IT^^,  can  oe  interpreted  as  the 

direct  and  the  perfect  ground  reflection  contributions  to  the  fields  for 

observation  points  in  region  1.  Therefore,  the  two  remaining  components 

„IL  , and  .IT,  , are  simply  the  correction  terms  to  the  perfect  ground 
0 hlx  0 hlz  v 3 

solution. 

2.4  Space-Domain  Representation  of  the  Vector  Potentials 


Since  for  most  practical  antenna  problems  the  observation  points  are 


located  above  ground,  in  this  section,  only  the  vector  potential  components 
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in  region  1 are  considered.  The  complete  transform  domain  expressions  for 
these  vector  potentials,  both  for  the  vertical  and  the  horizontal  current 
elements,  are  derived  in  the  previous  sections,  and  the  inverse  transform 
relation  (2.8b)  is  employed  to  obtain  the  corresponding  space-domain 
results . 

Bv  using  a spherical-type  change  of  variables,  viz., 

;'x  - r.,  sin  (?.,  cos  <p  - p cos 


/y  » r,  sin  0.,  sin  p.,  =■  p.,  sin 
^z  + h = r.,  cos  0 ^ * 2^ 


(2.35) 


oi  * -\  cos  z 

S * -\  s in  Z,  , ( 2 . 3b) 

the  incident  and  the  perfect  ground-reflection  components  of  the  vector 
potentials,  viz..  Equations  (2. 21a),(2. 21b)  and  (2 . 31a),(2.  31b)  can  be 
inverse  transformed  into  the  following  general  space-domain  form: 
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(2.37) 
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rr 


with  the  requirement  that  Imvk“  - \ 


0.  In  deriving  the  preceding 
equation,  the  following  identity  was  used 


cos(nX)J  (z) 
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(-1)  n ! .-.Izcos(t'-t) 

— . e cos 


2t 


cos  (nr ' ) dx  ' 


(2.38) 


th 


where  J is  the  n -order  Bessel  function.  Expression  C 2 . 5 7)  can  be 
n 

integrated  in  a closed  form  [15 1 to  vield: 


W * 1 1 exp  (-jk^r)/4rrr 


(2. 3°) 
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Equation  (2.39)  is  used  to  write  the  space-domain  expressions  for  the 
aforementioned  vector  potential  components  in  the  following  closed  forms: 


‘‘viz  vO 


Lv0  exp  (-jk1rl)/4irr1 


Viz  = ' rv0  exp  (_jkir2)/4rrr2 


nhlx  = Xh0  eXp  (-Jk1r1)/47Tr1 


- Ih0  exp  (-jk.  r )/4Ttr 


(2.40a) 


(2.40b) 


(2.41a) 


(2.41b) 


One  may  recognize  Equations  (2.40a)  and  (2.41a)  as  the  free-space 
Green's  function  solution  for  the  current  element  P^,  while  (2.40b)  and 
(2.41b)  are  the  corresponding  perfect  ground  contributions,  i.e.,  source 


at  image  point  P.,  , to  the  fields  at  observation  points  in  region  1. 

Inverting  the  remaining  vector  potential  components  in  region  1 and 
after  some  manipulations,  one  obtains: 


1 1 


popularly  known  .is  cho  Sommerfe  Id  integrals  | 1 | and  cannot  be 
expressed  in  a closed  form.  Efficient  numerical  evaluation  of  these 
integrals,  which  have  been  recently  verified  [?]  - [Is],  represents 
the  major  task  in  analysing  antenna  structures  over  lossy  ground. 

The  Sommer t eld  integral  can  take  several  forms.  For  this  work, 
the  forms  containing  Hankel  functions  in  their  integrands  are  preferred. 


since  the  integrals  appear  to  be  numerically  more  tractable.  Incor- 
porating the  well-known  identities  between  Bessel  and  Hankel  functions. 
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O'  hlz 


I.Ak  r _ / ' ] 

hO  1 I , 2 _ ,.  cos  t,  - /ic-sin 

— cos  sin  t,  cos  t. 

47T  _ j 


„ ( 2_ 
>s  E,  + K^-sin  t. 


/ -i \ -jk  z^cbst 

• (kjP,  sin  O e " d£ 


(2.49) 


where  Che  integration  path  F is  depicted  in  Figure  2,  on  which  the 
following  conditions  must  hold 

Im( cos  C)  < 0 (2.50a) 


/ 2 

Im(/<  - sin  5)  < 0 


(2.50b) 


In  the  following  chapter,  an  efficient  numerical  integration  scheme 
is  discussed  for  evaluating  the  Sommerfeld  integrals  via  the  steepest 
descent  path  integration  technique.  An  additional  approach  is  introduced 
in  Chapter  4 in  which  the  transform  domain  of  the  troublesome  vector 
potential  components  are  initially  approximated  such  that  the  Sommerfeld 
infinite  integrals  do  not  appear  in  their  space-domain  expressions.  These 
two  techniques  plus  the  asymptotic  solution  to  the  Sommerfeld  integrals 
are  compared,  and  several  numerical  results  are  included  to  demonstrate 
the  efficiency  and  the  accuracy  of  each  procedure. 
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3.  EXACT  AND  ASYMPTOTIC  EVALUATION  OF  THE  SOMMERFELD  INTEGRALS 

The  three  correction  vector  potential  components  derived  in  the 
previous  chapter,  viz..  Equations  (2 . 47)- ( 2 . 49)  , cannot  be  expressed  in 
a closed  form  and  contain  a certain  class  of  infinite  integrals  known  as 
the  Somnierfeld  integrals.  In  this  chapter,  an  efficient  technique  is 
presented  for  numerically  evaluating  the  aforementioned  integrals  by 
deforming  the  integration  contour  T into  the  steepest  descent  path  (SDP). 
The  asymptotic  approximations  to  these  Sommerfeld  integrals  are  also 
obtained  and  compared  with  the  numerical  integration  results.  Several 
numerical  examples  are  included  throughout  the  chapter,  which  demonstrate 
the  accuracy  and  the  efficiency  of  the  steepest  descent  integration 
technique  and  define  the  useful  range  of  the  asymptotic  expressions. 

3. 1 Steepest  Descent  Path  (SDP)  Integration 

It  can  be  readily  shown  chat  all  three  vector  potential  components 
given  in  (2.47)  - (2.49)  can  be  expressed  in  the  following  general  form: 

u = | P(S)  exp  jjjk^r,  cos  (£  - 6.,)"]  d£  (3.1) 

r 

where  P(£)  is  a relatively  slowly  varying  function  of  £.  The  point  £ = 0. 
is  a saddle  point  of  Equation  (3.1),  and  the  path  F can  be  deformed  to  a 
steepest  descent  path  (SDP)  by  enforcing  the  condition  Re [ cos (5  - 0->)  ] = 1. 
(At  this  point  it  is  assumed  that  no  poles  or  branch  points  of  P(%)  are 
intercepted.)  On  SDP,  the  following  convenient  change  of  variable  is 
introduced : 

cos  (£  - 0„)  = 1 - jt”  (3.2) 


i 


t 

i 
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where  t is  a real  parameter  ranging  from  -*»  to  +">.  Equation  (3.2)  can 


he  used  to  explicitly  define  the  steepest  descent  path  as: 

£gDp  “ t 77 1 2 + I Ln  (t"  + J + |t|  *t“  + 2])  I + d.,  : t * 0 , 


(3.3) 


or  if  one  needs  to  separate  the  real  and  the  imaginary  behaviors  of  the 


steepest  descent  path,  the  following  equivalent  form  can  be  derived 


( , r 2 nr~ " 

) - 1 - c + /t  + 4 . , .-It 

i ( cos  : + cosh 

\ L - L 


. -1  t“  + * t + 4 V . , 
sh  - / + d . 


t < 0 


(3.4) 


where  the  Lnverse  cosine  function  is  assumed  to  be  between  0 and  i/2 
and  the  inverse  hyperbolic  function  is  defined  as 


cosh  T - Ln  H'  + - 1 


(3.5) 


The  SDP  can  be  traced  in  the  f.-plane,  bv  using  Equations  (3.3’>  or  (3. 41 
as  a function  of  t and  0.,  (see  Figure  3).  Applying  the  change  of 
variable  in  Equation  (3.2)  Co  the  integral  expression  in  Equation  (3.1), 
we  obtain 


u - expi-jk^r.,  - Jn/4)/(2»2  n)  j Q(t)  exp  (~k^r,t“)  dt  (3.ti) 


where 


r « - e. 

Q(t)  * P(£)  sec  — 

L 


rr  1 1!  / 4 


J -e  2j)‘l/J  P(5)  it  • (3.7) 


Using  the  general  form  obtained  in  ( 3 . n ) , the  vector  potential  expressions 
in  f 2—4 7 ' - (2.41))  can  he  formulated  on  the  steepest  descent  path  as  a 


function  of  real  variable  t,  namelv. 
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where 


R^(t)  » (2n)  ^k  exp  (-Jk^r.,)  (t~  + 2j)  ^“[sin  t,  cos  t.  exp  (jk^p,  sin  .\)1 

• H|J\k.p0  sin  ^)  | ; 1-0,1  , (3.11) 

t i _ r 

Is  the  Hankel  function  of  the  i 1 order  and  of  second  kind,  and  i 
is  expressed  as  a function  of  t in  Equation  (3.3). 

The  relations  expressed  in  Equations  (3.8)  - (5.10)  are  exact  if  no 
poles  or  branch  points  are  Intercepted  under  the  path  deformation.  1'he 
detailed  discussions  of  the  possible  pole  and  branch-cut  cont r i but  ions  are 
presented  in  the  following  sections.  It  should  also  be  noted  that  the 
apparent  singularity  of  the  Hankel  functions  at  rV,  - 0 in  the  Sommerteld 
Integrals  is  overcome  bv  the  remaining  terms  in  the  integrands.  In  Appendix  1, 
equivalent  versions  of  Equations  (5.8)  - (3.10)  are  derived  tor  - 0 
for  numerical  integration  purposes. 

Bv  observing  the  integral  expressions  in  (5.8)  - (.5.105,  it  is 
apparent  that  for  large  k^r,  the  effective  integration  interval  will  be 
quite  small  and  contain  relatively  few  oscillations  (see  Figure  *1. 
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Fortunately,  at  Lower  frequencies  (smaller  k^r.,),  the  oscillatory  terms 
present  in  the  integrands  are  scaled  accordingly  and, even  though  the 
effective  integration  interval  increases,  the  number  of  oscillations  in 
the  integrand  does  not  increase  appreciably  (Figure  5).  This  fact  allows 
one  to  integrate  Equations  (3.8)  - (3.10)  numerically  by  employing  an 
efficient  Gaussian  quadrature  integration  routine  for  a wide  range  of 
parameters.  Later  in  this  chapter,  the  efficiency  and  the  accuracy  of 
the  above  numerical  integration  procedure  are  demonstrated  and  compared 
to  the  other  available  techniques. 

3 . 2 Branch-Cut  Contribution 

The  expressions  derived  in  the  previous  section  are  valid  only  when 
no  singularities  are  intercepted  during  the  steepest  descent  path  deformation. 
In  order  to  locate  the  poles  and  the  branch  points,  one  must  consider  the 
following  physical  constraints: 

a)  0 <_  0^  < ir/2,  since  Equations  (3.8)  - (3.10)  are  valid  only  for 
observation  points  above  ground; 

b)  Re(K)  >_  1 and  Im(ic)  <_  0,  since  < = -ja/(u)£Q); 

c)  - -y  < Re(5)  < it  on  the  SDP  (see  Figure  3) . 

Furthermore,  since  cos  (5)  is  a single-valued  function,  condition  (2.50a) 
does  not  have  to  be  satisfied  during  the  path  deformation.  Condition 
(2.50b)  is  used  to  define  an  upper-  and  a lower-Riemann  sheet  in  the 
5-plane  in  which  this  condition  is  satisfied  in  the  upper  and  violated  in 
the  lower  sheet. 

Equations  (3.8)  - (3.10)  have  the  same  branch  points  satisfying 


< - sin 


(3. 12) 


t »■  HIM 


Low-t requeney  examples  at  c lie 
h ■ im,  : * 10,  and  j « .01  mhos  m. 


Figure  3 


For  these  oases 


If  one  considers  the  physical  constraints  discussed  earlier,  only  the 
following  two  branch  point  solutions  are  of  importance  (see  Figure  o)  : 


n/2  i j Ln  (ik  + »<  - 1) 


(3.13) 


The  corresponding  branch  cuts  of  Equation  t.3.13)  which  satisfy  the  relation 


Ira(»,<-sin~S)“0  , 


(3.14) 


as  depicted  in  Figure  6,  are  the  boundaries  through  which  the  integration 
path  will  travel  to  and  from  the  two  Riemann  sheets  defined  earlier. 

Since  the  steepest  descent  path  defined  in  Equation  (3.3)  is  independent 
of  < , one  can  easily  demonstrate  that  for  0^0,  < dO"  only  the  branch 
point  with  the  upper  sign  can  be  captured  bv  the  SDP  deformation  (Figure  7). 
Therefore,  one  can  allow  the  SDP  to  enter  the  lower  sheet  onlv  when  the 
lower  branch  cut,  corresponding  to  the  lower  sign  of  Equation  t.3.13),  is 
intercepted,  since  the  path  will  always  intercept  the  lower  cut  at  an 
additional  point  forcing  it  to  return  to  the  upper  Riemann  sheet  (see 
Figure  7).  A branch-cut  integration,  however,  is  performed  around  the  upper 
branch  cut  whenever  it  is  intercepted  in  order  to  remain  the  proper  sheet. 
The  branch  cut  in  the  upper-half  plane  as  a function  of  a positive  real 
parameter  3 can  be  expressed  as 


it/2  + j Ln  (/<  -3  + »'K  — 1 - ) 


(3. 1 5a) 


or,  equivalent  lvt 


(3.  15b) 


-j  /<  -1-3* 


(3.  1 5c) 


Bv  applying  the  change  of  variable  in  Equation  (5.  I5')  to  the  Sommer f eld 


Branch  Point  Loci 


Brancli~poi.nl  and  branch-cut  loci  as  a function  of  ground  parameters  L and  o/f 


Jb 


integrals  expressed  in  tE.s73  - (E.sO)  and  integrating  around  the  branch 
cut,  the  following  contributions  are  obtained: 
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13.183 

where  3 ” 8^  is  the  crossing  point  of  the  branch  cut  and  the  SDP  in  the 
'.-place.  By  expanding  Equation  (3.J3  and  by  employing  the  (3.153 
relations,  the  following  conditions  are  obtained  for  8, 


sin  0,  Re (A 3 + cos  8,  tmlB3  * 1 


sin  8,  lm(A)  + cos  8,  Re(B3  1 » t. 


where  the  complex  numbers  A and  B are  defined  as 


A » » < - 3~  ; Re(A3  ' 0 , ImiA3  ^ 0 


8 = > < - 1-3“  ; Re(83  ^ 0 . lm(B3  0 


t 3.  10.,  3 


(3.  1 9b 3 


(3 . 80a 3 


1 3 . 30b 3 


In  (3.193,  t * t^  defines  the  point  at  which  the  SDP,  corresponding  to 
the  observation  angle  5,.  intercepts  the  branch  cut  of  Equation  y 8 . 1^3  .,t 


After  some  algebraic  manipulations.  Equation  (3.19a)  can  be 


further  simplified  to 

8 , “ sin  ' ilRe'iA)  + lm'lB)]  ' ~)  - tan  ^ (pr^TT") 
_ Ke 1 A ) 


13.21) 


from  which  8 . . the  observation  angle  at  which  SDP  will  pass  through  the 
min 

branch  point,  can  be  computed  bv  setting  6,  * 0.  Since  8 , is  onlv  a 

1 min 

function  of  c,  Figure  8 is  constructed  to  show  its  variations  as  a function 
of  the  ground  parameters  and  the  frequency.  Therefore,  whenever  the 


observation  angle  8 , satisfies  the  condition 


13.22) 


the  branch-cut  contributions  in  l3.1t>)  - i3.18)  are  to  be  added  to  their 
respective  SDP  vector  potential  formulations  expressed  in  i3.8)  - 13.10). 

It  should  be  pointed  out  that  once  the  condition  13.22)  is  met,  the 
branch-cut  integration  limit  d ^ can  be  computed  numerically  bv  iterating  on 
Equation  id. 21).  Also,  because  of  the  branch-cut  intercept  ion , the  SDP 
integrand  will  be  discontinuous  at  point  t = t^,  which  is  readily  computed 
bv  substituting  the  value  of  8 into  i3.1db). 

Fortunately,  in  many  cases,  the  branch-cut  contributions  are  several 
orders  of  magnitude  smaller  than  the  SOP  integral  value  and  can  be 
ignored  [In],  Therefore,  it  is  necessary  to  introduce  a condition  for 
which  one  can  Ignore  the  branch-cut  integration  and  thereby  compute  the 
vector  potentials  more  efficiently.  Phis  task  can  be  accomplished  by 
initially  considering  the  exp  l-k^r.t')  term  present  in  all  three  vector 
'otential  integrands  shown  in  i3.S)  - id.lO''.  If  no  poles  are  present 
•n  >r  near  the  contour,  a finite  integration  in  the  interval 
a 1 ' 

ir— ->  ' 13.23) 

— tv , r , 


13.23) 


will  result  in  an  error  on  the  order  of  0.012  as  compared  to  the  full 
infinite  integration.  By  examining  the  branch-cut  loci  (Figure  6)  and 
the  SDP  behavior  (Figure  3)  in  the  £-plane,  an  assumption  can  be  made 
which  states  that  if  the  branch  cut  intercepts  the  SDP  inside  the  finite 
integration  interval  defined  in  (3.23),  then  the  branch-cut  contribution 
is  not  negligible.  Or  equivalently,  the  branch-cut  contribution  requires 
an  additional  condition,  namely. 


t,  < t 
1 — max 


3(1^)' 


(3.24) 


where  t = t^  is  the  SDP  and  branch-cut  intercept  defined  earlier. 

In  summary,  condition  (3.22)  signals  the  capture  of  the  branch 
point  during  the  path  deformation.  If  captured,  condition  (3.24)  is 


used  to  decide  whether  the  branch-cut  integration  can  be  ignored  or  not. 


Vv 

\i 

V 


Table  3.1  is  constructed  to  verify  the  validity  of  the  assumptions  which 
led  to  Equation  (3.24)  by  comparing  the  branch-cut  values  with  the  SDP 


integration  results  for  a wide  range  of  parameters. 


(■ 

3.3  Pole  Contribution 


Unlike  the  branch  points 


which  exist  in  all  three  of  the  correction 


vector  potential  components,  the  poles  only  exist  in  the  vertical 


components  of  the  vector  potentials,  viz 


the  previous  section,  only  the  following  two  poles  need  to  be  considered 


Fortunately,  one  can  verify  that  for  0 y o/f  •-  M,  0 y 0,  » 90“ , and 
by  considering  the  correct  Rlemann  sheets  (see  Figure  7),  the  poles  In 
liquation  ().2b)  will  not  be  captured  bv  the  SDP  and  the  residue  contribution 
is  not  needed.  However,  under  extreme  circumstances,  I.e.,  0.,  90°  and 

small  i t,  a pole  can  come  close  to  the  SDP.  In  this  case,  a higher-order 
Gaussian  quadrature  Integration  routine  Is  required  and  possibly  the 
effective  Integration  Interval  In  (3.21)  should  be  expanded. 

1 . 4 Asvmptot 1c  Approx (mat Ion 

Che  steepest  descent  formulation  of  the  previous  sections  can  be  used 
to  derive  an  asymptotic  expansion  for  the  Sommert'eld  integrals  in  terms 
of  Inverse  powers  of  k^r,.  In  Appendix  II,  a general  discussion  Is 
presented  for  asvmptot  lea l lv  evaluating  the  integral  in  (3.1),  and  as  an 
example,  the  first  two  asymptotic  expansion  terms  lor  the  vector  potential 
expressions  In  (2.47)  - (2.49)  are  derived.  The  first  terms,  better  known 
as  the  Fresnel's  reflection  coefficient  method  (ROl)  approximations,  are 
shown  to  be: 
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The  above  RCM  expressions  have  been  extensively  used  for  high-frequency 
antenna  applications  large)  with  surprisingly  accurate  results, 

e.g.,  [13],  [14].  This  success  has  been  mainly  due  to  the  fact  that 
for  observation  points  away  from  the  interface,  the  remaining  vector 
potential  components  present  in  (2.20)  and  (2.30)  tend  to  dominate  over 
the  RCM  components  thereby  reducing  the  net  error  in  the  total  vector 
potential  values  computed.  Therefore,  little  is  gained  by  adding  the 
complicated  second  terms  in  the  asymptotic  expansion  of  the  vector 
potentials,  derived  in  Appendix  II,  since  the  accuracy  of  the  total 
vector  potential  will  not  be  affected  appreciably. 

An  error  of  5%  or  less  can  be  expected  in  the  RCM  vector  potential 
components  shown  in  (3.27)  - (3.29)  when 

> 10  , (3.30) 

and  as  long  as  the  branch-cut  conditions  given  in  the  previous  section 
are  not  violated  (see  Tables  3.2  - 3.4).  These  conditions  are  more 
general  than  the  one  proposed  by  Sarkar  [ 7 ] in  defining  the  useful 
range  of  the  RCM  expressions.  The  second  terms  in  the  asymptotic  expansion, 
as  demonstrated  in  Tables  3.2  - 3.4,  only  offer  a slight  improvement  in 
the  accuracy  of  the  JI  expression  in  the  region  where  the  branch- 

cut  contribution  is  negligible.  The  remaining  vector  potential  components, 
however,  do  not  benefit  from  the  2nc*  term  in  the  asymptotic  expansions, 
possibly  because  they  contain  a pole  in  their  Sommerfeld  integrals. 

Figures  10  - 12  also  compare  the  RCM,  two-term  asymptotic  expansion  and 
the  exact  integration  values  of  the  correction  vector  potential  components 
for  a typical  half-space  as  a function  of  k^r,,.  Finally,  in  order  to 
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demonstrate  the  branch-cut  effect.  Figure  13  is  included  to  show  the 
minimum  k^r,  contours  in  the  ic-plane  above  which  one  can  neglect  the 
possible  branch-cut  contribution  (based  on  condition  (3.34)). 

The  computing  time,  of  course,  is  the  main  reason  why  the  RCM 
approximut ions  are  used  whenever  possible.  For  example,  the  execution 
time  on  a Cyber  175  computer  for  computing  the  three  correction  vector 
potentials  is  » 1 msec  for  the  RCM  approximation  and  as  high  as  50-100  msec 
for  the  SDP  integration  technique  presented  in  this  chapter.  It  should 
be  pointed  out,  however,  that  the  SDP  integration  is  an  order  of  magnitude 
faster  than  many  of  the  previously  reported  integration  techniques 
13]  - [9],  [12],  making  it  more  suitable  for  a much  wider  range  of  low- 
frequencv  applications. 
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4.  COMPUTATIONS  OF  THE  VECTOR  POTENTIALS  WITHOUT 
SOMMERFELD  INTEGRATION 


The  conventional  approach  for  analyzing  antenna  problems  over  lossy 
ground  involves  repeated  evaluations  of  the  Somraerfeld  integrals  appearing 
in  the  expressions  for  the  vector  potentials,  viz..  Equations  (3.8)  - (3.10). 
Even  though  the  SDP  procedure  discussed  in  the  previous  chapter  is  an 
efficient  integration  technique,  the  computation  time  needed  for  evaluating 
these  integrals  severely  limits  the  physical  dimensions  of  the  antenna  problems 
being  analyzed.  In  this  chapter,  a novel  approach  is  introduced  by  initially 
approximating  the  well-behaved  transform  domain  expressions  of  these  integrals, 
derived  in  Chapter  2,  and  then  performing  the  inverse  transform  operations 
via  a set  of  exact  identities.  The  resulting  space  domain  expressions  are 
valid  for  a wide  range  of  parameters,  with  their  computation  times  being 
comparable  to  those  of  the  RCM  approximation. 

4 . 1 Transform  Domain  Expressions 

The  transform  domain  expressions  for  the  three  correction  vector  potential 
components  shown  in  Equations  (2.21c),  (2.31c),  and  (2.32)  can  be  rewritten 
in  the  following  forms: 
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where  z^  = z + h,  and  y 


h 
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defined  in  (2.15),  are  explicitly  expressed  as 


Itn (y^  ) <_  0 
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(tck,  - ot“  - 3 ) ~ 


; IraCy.,)  < 0 


In  addition,  the  free-space  Creen’s  function  g can  also  be  expressed  in 


(4.4a) 


(4.4b) 


the  transform  domain  (see  Equations  (2.41b)  and  (2.31b))  yielding  the 
following  transform  pair: 


g(0t,  £,  Z,)  = exp  (-jY12,) 


(4.5a) 


g(x,  y,  z0)  = exp  (-jk  r ,)/4iTr.,  ; r. 


2 2 7 1 / 7 

[x“  + y + z ] /_ 


(4.5b) 


By  applying  a successive  3/3z  operator  to  (4.5),  an  infinite  set  of  trans- 
form pairs  is  obtained,  viz. , 


Qn  a y"  L exp  (-JY1z2) 


(4.6a) 


?mn+1  JL 


g(x,  y,  z ) 


n = 0,1,2,. 


(4.6b) 


where  Q is  expressed  in  a closed  form  for  all  n and  can  be  numerically 

evaluated  quite  rapidly  (see  Appendix  III). 

An  examination  of  the  Fourier  transformed  vector  potentials  in  (4.1)  - 

(4.3)  and  the  in  (4.6a)  reveals  two  important  and  useful  properties. 

First,  all  of  the  expressions  have  an  identical  z-variation  term  that 

corresponds  to  a space-domain  solution  emanating  at  the  image  point  P,. 

Second,  it  is  apparent  that  all  of  the  equations  are  well-behaved  in  the 

Fourier  domain  and  decay  exponentially  to  zero  outside  the  circle 
2 2 2 

+ 3“  k~.  These  properties  give  rise  to  the  possibility  of  performing 

the  inverse  transtorm  operation  on  Equations  14.1)  - r. 4 . 3 ) via  the  use  of 
the  (4.6)  identities.  The  major  obstacle  to  such  a procedure  is  the  existent 
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of  y,,  in  these  expressions  which  is  overcome  in  the  following  section  by 
an  appropriate  approximation. 

4.2  Approximating  y.,  and  Space-Domain  Results 


All  three  Fourier  transform  domain  expressions  in  (4.1)  - (4.3)  can 
be  put  into  the  following  general  form 

(4.7) 


0n  - f (yl » yn)  exp  (-jy^z7) 


As  mentioned  in  the  previous  section,  an  important  property  of  Equation  (4.7) 

is  the  exponential  term  which  rapidly  decays  to  zero  outside  the  circle 
2 2 2 

a”  + 8“  = k“  (also  see  Figure  14).  This  fact  enables  one  to  replace  y, , 
defined  in  (4.4b),  by  the  first  term  of  its  Taylor's  series  expansion  y^, 
that  is. 
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(4.8a) 

(4:8b) 


For  most  practical  values  of  v,  i.e.,  !k|  > 10,  the  approximation  in  (4.8b) 


is  excellent  inside  the  circle  oT  + 3 k“.  Fortunately,  as  demonstrated 


in  Figure  14,  the  decaying  exponential  in  (4.7)  can  easily  overcome  the 

1 1 9 


errors  introduced  by  (4.3b)  in  the  region  + 8 k“  thereby  making  this 


a valid  approximation  throughout  the  cx8-plane  for  a wide  range  of  k and 
z^  parameters,  namely. 


(4.9) 


0n  = f(y1>  y2)  exp  (-jkLz2) 

where  the  bar  on  top  represents  the  approximate  quantities.  It  should  be 
noted  that  Kuo  and  Met  [12,  Eq.  8[  have  recently  verified  and  used  the 
aforementioned  approximation  for  simplifying  the  space  domain  expressions 


: 


Mil 


] 

. ...J 


while  Chang  et  al.  [b]  have  employed  an  equivalent  approximation  for 
representing  the  vector  potential  in  the  vertical  current  element  problem 
in  terms  of  an  incomplete  Hankel  function.  The  approximate  transform  domain 
vector  potential  quantities  can  now  be  treated  individuall'’  in  order  to 
obtain  their  respective  transform  domain  expressions. 

•4.1.1  Approximation  for  the  horizontal  component  Applying 

the  approximation  introduced  in  Equation  (4.8b)  to  the  horizontal  correction 
vector  potential  component  in  (4.1),  one  obtains 
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(4. 10) 


and  direct  application  of  the  (4.b)  Identities  results  in  the  following 
closed-form  space  domain  expression: 

-II. 
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(.4.11) 


Note  that  the  various  partials  of  g are  explicitly  derived  in  Appendix  III. 
The  useful  range  of  and  the  comparison  of  Its  accuracy  and  efficiency 

with  other  available  techniques  are  discussed  in  the  next  section. 

4.1.1  Approximation  for  the  vertical  components  (^!I  ^ and 

Again  the  variable  y,  can  be  eliminated  in  (4.1)  and  (4.3)  by  applying  the 
approximation  in  (4.8b),  that  is. 
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where  the  constant  c is  defined  as 


k . / » V 


(4. 14) 


Application  of  the  identities  in  Equation  (4.6)  results  in  the  following 
space  domain  forms: 
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a unique  solution  to  the  (4.18)  differential  equation.  Since  no  simple 
starting  point  can  be  found  at  which  S(p,  z\)  is  known,  it  is  assumed  that 
is  large  enough  so  that  the  RCM  approximation  in  Equation  (3.17)  is 


applicable  (see  Figure  15),  that  is. 
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(4.20) 


where  angle  0,,  is  defined  in  Figure  15.  A simiLar  procedure  can  be  used 

8 - 

for  computing  ^ S,  required  for  Qn  in  Equation  (4.15),  by  simply 
differentiating  (4.19)  with  respect  to  x: 
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(4.21) 


and  by  using  the  asymptotic  expression  in  (3.29),  the  following  initial 
value  at  point  z^  is  obtained: 


it  S(p’  Z2> 
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In  summary,  the  procedure  for  computing  the  vertical  vector  potential 
components  at  the  point  (p,  z,,),  as  demonstrated  in  Figure  15,  is  to  start 
from  a higher  observation  point  (p,  z',).  at  which  the  RCM  expressions  are 
valid,  and  simply  use  Equation  (4.19)  or  (4.21)  to  integrate  down  to 
(p,  z,,).  This  procedure  effectivelv  computes  a correction  term  for  the 
RCM  approximations  of  the  vertical  vector  potential  components  in  the 
region  where  RCM  alone  breaks  down.  The  useful  range  and  the  accuracv 
of  this  procedure  are  discussed  in  the  following  section. 


4. 3 Error  Estimation  tor  the  Approximate  Expressions 

In  the  previous  section,  only  the  approximation  in  (4.8)  was  made  to 
derive  the  horizontal  vector  potential  component  in  Equation  (4.11).  The 
approximated  term  T in  this  vector  potential  component,  for  convenience, 
can  be  expressed  in  the  following  forms: 


T = Y2  exp 

(_jYlZ2) 

(4.23a) 

T = y2  exp 

(_JY1Z2) 

(4.23b) 

T = -Ik^*: 

JI  8tp’  z2) 

(4.23c) 

The  error  introduced  by  (4.23b)  in  the  transform  domain  is  simplv 

AT  = T - T = (y2  - y2)  exp  (-jY^,)  . (4.24) 

and  again  referring  to  the  Taylor's  series  expansion  for  y.  in  (4.8a), 

2 2 2,  i 

which  is  convergent  for  a +8  < k^;<!,  one  can  assume  AT  to  be 

proportional  to  the  second  term  of  the  expansion,  namely, 

2 2 

AT  = -<a"  -+  3 ) exp  (-JY,*,)  • (4.25) 

2k./<  - 

At  the  same  time,  a constraint  is  needed  for  in  order  to  enforce  the 
exponential  term  in  (4.25)  to  decay  sufficiently,  i.e.,  to  1*,  when 
ot“  + 3“  = k“|<|.  This  important  condition,  which  allows  one  to  use  (4.25) 
for  error  analysis,  can  be  expressed  in  the  following  simple  form: 

k z,  > ---J  . (4.26) 
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Bv  substituting  \ + 3“  from  (4.4a)  into  Equation  (4.25),  one  can  obtain 

AT  directly  via  the  (4.6)  identities,  namely, 

kl  3 1 

)T  = — Yi  f*(P»  + 2:  — 3 8(0,  z,)  . (4.27) 
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The  space  domain  expressions  in  (4.23c)  and  (4.27)  can  be  used 
R^,  the  relative  error  in  the  magnitude  of  T,  as 


AT'  < 
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and  by  using  the  explicit  expansions  provided  in  Appendix  III, 
relation  can  be  further  reduced  to 


to  define 


(4.28) 

the  above 


Equation  (4.29)  is  a useful  upper  bound  for  the  error  existing  in  the  computed 
magnitude  of  gll^  . The  first  term  in  the  aforementioned  equation  sets  a 
limit  on  the  minimum  required  [k|  value;  for  example,  the  error  condition 
R^,  < 10%  requires  |k|  > 10,  while  R^  < 5%  will  require  |<|  > 20.  Figures  16 
and  17  demonstrate  the  valid  regions  in  the  (r„,  0.)  plane  as  a function  of 
|<|  for  two  useful  error  conditions,  namely,  R^,  = 5%  and  10%,  respectively. 
Table  4.1  is  included  to  numericallv  verify  the  conditions  (4.26)  and  (4.29) 
by  comparing  the  approximate  horizontal  vector  potential  values  in  (4.11) 
with  the  corresponding  exact  SDP  integration  results  of  Chapter  3.  As 
clearly  demonstrated  in  this  table,  the  actual  computed  errors  in  the 
magnitude  are  consistently  below  the  R values  obtained  from  Equation  (4.29), 
for  the  various  (k,  r0,  0,)  combinations  tested. 

Unfortunately,  because  of  the  complicated  nature  of  the  approximate 
transform  expressions  (4.12)  - (4.13),  one  cannot  derive  a simple  error 
condition  for  the  vertical  potential  components,  as  was  derived  for  the 
horizontal  case  (Equation  (4.29)).  However,  because  of  the  similar  nature 


53 


of  the  approximation  for  all  three  of  the  correction  vector  potential 
components,  one  can  assume  that  conditions  14.26)  and  (4.29)  also  apply 


to  the  _ and  q^v^z  expressions  shown  in  Equations  (4.15)  - (4.1b). 

This  assumption  is  verified  bv  the  results  in  Tables  4.2  and  4.3  in 
which  the  vertical-correction- vector  potential  components  are  computed 
by  using  k^r',  = 10  as  a starting  point.  Even  though  these  ROM  initial 
values,  in  general,  contain  about  51  error,  the  net  errors  in  all  of  the 
cases  are  well  below  the  predicted  ones  (Equation  (4.29)).  Table  4.4  is 
also  included  to  demonstrate  the  stability  of  the  final  vertical  vector 
potential  value  as  a function  of  starting  point  z,'.  Note  that  when 
k^r',  > 10,  the  error  in  the  final  results  stays  within  acceptable  levels 
of  tiie  exact  values  primarily  because  the  RCM  approximation  used  to  compute 
the  starting  point  becomes  reasonably  accurate  in  this  region. 

Studying  the  results  of  this  section,  one  can  conclude  that  the 
proposed  approximate  formulas,  i.e.,  (.4.11),  (4.13)  and  (4. In),  can  be 


employed  in  most  practical  antenna  problems  as  long  as  V > 10  and 
condition  (4.26)  is  satisfied. 

4 . ->  Advantages  of  the  Approximate  Expressions 

The  major  advantage  of  the  approximate  technique  introduced  in  the 
chapter  is,  of  course,  the  computation  of  the  vector  potentials.  For 
example,  on  the  Cyber-175  computer,  the  computation  time  for  computing 
the  three  vector  potential  expressions  is  ~ l msec  for  the  RCM  technique, 
- 5 msec  for  the  approximate  technique  presented  in  this  chapter,  and 
finally  - 50-100  msec  for  the  SDP  integration  technique  discussed  in 
Chapter  3.  Note  that  even  though  the  approximate  technique  is  slightly 
more  time-consuming  than  the  RCM  method,  it  can  be  applied  to  a much 


Table  4.1  Comparison  of  the  KCM,  exact,  and  the  approximate  evaluations  of  II  . . For  this  example,  f - 10  MHz, 
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wider  range  of  parameters  making  it  suitable  for  most  practical  antenna 
problems. 

* Another  advantage  of  this  technique  is  that  the  horizontal  vector 
potential  component  is  expressed  in  a simple  closed  form  (Equation  4.11), 

while  Kuo  and  Mei  [12]  have  obtained  yet  another  infinite  integral  form 
by  applying  virtually  the  same  approximation  as  the  one  in  Equation  (4.8b) 
to  the  Sommerfeld  integrals.  Also  the  finite  integrations  needed  for 
the  vertical  vector  potential  approximations,  viz..  Equations  (4.19)  and 
(4.21),  will  never  have  a singularity  in  their  integration  interval  thus 
making  the  integrands  well-behaved.  In  addition,  these  integrands  are 
independent  of  z7,  thereby  making  it  possible  to  compute  the  vertical 
vector  potential  values  along  the  vertical  interval  z7z7  by  a single 
integration,  i.e.,  by  using  the  newly  computed  value  of  the  vector 
potential  as  an  initial  value  for  computing  the  next  point  on  the 
interval.  This  procedure,  when  needed,  can  appreciably  improve  the 
overall  efficiency  of  the  technique. 

Finally,  since  the  evaluation  of  the  E-  and  the  H-fields  is  of  major 
importance,  the  approximate  technique  presented  in  this  chapter  also 
provides  a computationally  efficient  formulation  for  computing  the  various 
field  components.  This  can  be  simply  demonstrated  by  observing  the 
matrix  Equations  (2.9)  and  (2.10).  The  various  mixed  partial  derivatives 
of  and  II  , needed  for  the  field  computations,  can  be  easily  expressed 
in  terms  of  the  various  mixed  partial  derivatives  of  g(p,  z7) , which  in 
turn  are  expressed  in  a computationally  efficient  closed  form  in  Appendix  III. 

The  following  chapter  demonstrates  the  ability  of  these  approximate 
expressions  to  efficiently  analyze  several  antenna  structures  radiating 
over  a lossy  half-space. 
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5.  WIRE  ANTENNAS  RADIATING  OVER  A LOSSY  HALF-SPACE 

The  approximate  field  solution  to  the  current  element  problem  radiating 
over  a lossy  half-space,  developed  in  the  previous  chapter,  in  conjunction 
with  the  method  of  moments  [17]  can  be  employed  to  analyze  a wide  variety 
of  thin-wire  antenna  problems  radiating  over  a lossy  half-space.  Initially, 
in  this  chapter,  a general  integral  equation  is  derived  containing  the  unknown 
antenna  current.  The  method  of  moments  is  then  used  to  reduce  the  integral 
equation  into  a numerically  manageable  matrix  form,  and  finally,  a digital 
computer  program  is  developed  for  computing  the  antenna  currents,  impedance, 
and  far-field  patterns,  given  a specified  antenna  geometry.  A number  of 
simple  antenna  structures  are  considered  and  their  behaviours  are  numerically 
predicted  and  presented  in  the  final  sections. 

5 . 1 Antenna  Integral  Equation 

Figure  18  depicts  the  geometry  of  an  arbitrary  wire  antenna  over  a 
lossy  half-space,  with  (r  ,9  ,j>  ) defining  a point  on  the  antenna-axis. 

3 3 3 

For  simplicity,  it  is  assumed  that  the  antenna  is  entirely  in  the  xz-plane 
(j>  = 0),  since  the  field  computation  due  to  currents  in  the  v-direction  can 

3 

easily  be  handled  by  a digital  computer  program  by  simply  rotating  the 

xy-plane  by  90°  about  the  z-axis.  Assuming  that  the  antenna  is  excited 

+exc  ■> 

by  the  field  E (r  ),  and  having  a loading  function  A(r  ) ohms/meter, 

3 3 

one  can  write  a general  integral  equation  enforcing  the  total  E-field  a Lone 
the  antenna  equal  to  that  induced  by  the  possible  loading  function,  that 
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Ur)  ■ P“(t)  + i(r.)  ■ • t(7I>  + ■ U^K 


antenna 


= A(r  )I(r  ) , 

a a 


where  I(r  ) is  a unit  vector  along  the  antenna  direction  and  the  kernels 

ci 

G (r  ,r')  and  G,  (r  ,r')  are  the  E-fields  induced  at  point  r due  to  a one- 

V SL  2L  tl  SL  Si  cl 

ampere  electric  current  element  located  at  r^  and  oriented  in  the  z-  and 
the  x-directions , respectively.  The  matrix  Equation  (2.9),  which  formulates 
the  electric-field  components  in  terms  of  the  vector  potentials,  is  used 
to  write: 


2 = t>  n + d n , 

xx  z z 


where  the  vector  operators  D and  D are  defined  as: 

X 2 

D . i[k?  + -4  I + yt4-l  + 

x [1  ..  2 ^(3x3yj  [3x3z 


(5. 3a) 


D = x — — r — 
z 3x3z 


+ + 2k?  + -SL 

W3zJ  I 1 3z2 


(5.3b) 


By  using  the  general  formula  in  (5.2),  one  can  directly  write  the  two  kernels, 

G and  G,  , in  terms  of  the  incident,  perfect  reflection,  and  the  correction 
v h r 

vector  potential  components  (see  Equations  (2.20),  (2.30  . and  (2.32)),  namely, 


G (r  ,r')  = G*(r  ,r')  + Gr(r  ,r')  + G (r  ,r'), 
v a a v a a v a a Ov  a a 


where 


(5.4) 


Gv(ra’ea>  = (ju)£0)  Dz8(ra  " 


(5. 5a) 


Gr(r  , r 1 ) = ( juic_)  g(r  - r’  - 2hz) 
vaa  0 z a a 


(5. 5b) 


0Gv(ra’ra) 


D 1 (r  - r'  - 
z 0 viz  a a 


(5.5c) 
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and 


(1  (r  ,r')  = G*(r  ,r')  + oF(r  , r ' ) + .G,  (r  ,r')  , 

h a a h a a h a a 0 h a a 


(5.6) 


where 


= ^V'^x^a  - 


Bh(Ba’Ba}  = (j“e0)_lBx8(Ba  " Pa  ' 2h£) 


oV'a'^ 


[Bx  0* 


, . (r  - r1  - 2hz)  + D II,  , (r  - r1  - 2hz) 
hlx  a a zOhlza  a 


(5.7a) 

(5.7b) 


(5.7c) 


Th  - 1 


Note  that  in  the  above  equations,  h = r’  cos  (0  ) is  the  height  of  the  current 
source  above  the  half-plane  interface  and  g is  the  free-space  Green's  function 
defined  in  (4.5b). 

The  correction  vector  potential  formulation  of  Chapter  4,  namely. 

Equations  (4.11),  (4.15),  and  (4.16),  along  with  the  expansions  presented 

in  Appendix  III,  enable  one  to  compute  the  scattered  components  of  the  two 

kernels,  viz.,  G ,rtG  ,G,  , and  -G,,  without  difficulty.  However,  because  of 
v 0 v*  h’  Oh* 

the  singular  nature  of  g,  the  free-space  solution  of  the  kernels,  G^  and  G^, 
should  not  be  computed  directly.  Instead,  as  has  been  successfully  reported 
[11],  [14],  [18],  and  [19],  the  thin-wire  approximation  is  used  to  shift 
the  observation  point  r^  to  the  antenna  surface,  and  in  order  to  further 
smooth  out  the  singularities,  the  finite  difference  scheme  is  employed  to 
perform  the  and  operations  defined  in  (5.3). 

5. 2 Method  of  Moments 

As  developed  by  Harrington  [17],  the  method  of  moments  is  a convenient 
approximation  for  transforming  the  antenna  integral  equation  into  a numerically 
manageable  matrix  form.  In  this  work,  pulse-basis  and  delta-matching  functions 
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are  chosen  since  they  eliminate  the  need  for  integrating  the  kernels  £ 


and  G,  . The  number  of  unknown  patches  on  the  antenna  N should  be  large 
h 

enough  so  that  the  patch  length  d is  at  most  1/6  of  the  wavelength.  The 


approximated  current  along  the  antenna  is  therefore  represented  as: 

N 

i “ -ill  . (5.s: 

n»l  n n 

for  which  I is  an  unknown  constant  value  over  the  n1*1  patch  and  zero 
n r 

outside  of  it,  aiso  I is  a known  unit  vector  tangent  tc  the  3ntenna  at 
n 

the  center  of  the  n patch  (see  Figure  18).  Substituting  (5.8)  into  i 5 . 1 ) 
and  letting  subscripts  n » 1,2,3,***  denote  "evaluation  at  the  center 
of  the  n^1  patch,"  and  letting  [ f]  and  [Eexc]  be  column  vectors  containing 
the  current  and  the  tangential  excitation  field  values  at  successive 
patches,  one  finally  arrives  at 

[Eexc]  . -[Zimpim  + [A]  [II  , (5.9: 


where  (A)  is  a diagonal  matrix  with  elements  \^,Aof • • • ,A  ; and  [Z  I 

an  n x n square  matrix  with  its  it*1  row  and  j11*'  column  element  defined  as: 


Z = d 1 (r  ) • [z  • I(r  ,)G  (r  r ) + x • I(r  )G  (r  .,r  .)]  . (.5.10) 

ij  ai  aj  v ai  aj  aj  h ai  aj  1 

Note  that  i and  j also  refer  to  patch  numbers,  and  a is  the  patch  length. 

The  matrix  Equation  (5.9)  can  be  solved  for  the  unknown  currents  [I],  and  by 

replacing  the  excitation  E-field  in  terms  of  the  excitation  voltage  [V],  one 

arrives  at 


[l]  - [Yantl[V] 


(5. 11a) 


[YJnt)  . A_1i - [Zlmpl  + [a])"1  . 


(5.11b) 


Once  [Y  1 is  constructed  for  a given  structure,  the  antenna  currents  e, 
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be  directly  computed  for  a given  voltage  source  excitation.  As  was  discussed 

in  the  previous  section,  the  thin-wire  approximation  and  the  finite-dif ferenct 

scheme  are  empioved  for  evaluating  the  free-space  solution  components  of 

the  kerneis,  viz.,  G*  and  . A power  series,  derived  bv  Harrington  [20], 
v h 

is  used  for  the  thin-wire  approximation  computations , and  based  on  the 
conclusions  made  in  [11],  [14],  [18],  and  [19],  the  finite  difference 
parameter  6 is  chosen  to  be  equal  to  A/2,  e.g., 


3x~ 


f(xQ  + 5)  + f(xQ  - S)  - 2f(xQ) 


A/2  . 


(5.12) 


The  free-space  solution  obtained  by  using  the  above  approach  has  been  thor- 
oughly tested  and,  as  an  example,  the  generated  impedance  curves  shown  in 
Figures  19  and  20  agree  well  with  the  ones  reported  by  Jordan  et  al.  [21]. 

5 . J Far-Fleld  Radiation  Pattern 

The  RCM  expressions,  shown  in  Equations  (3.27)-  (3.29),  are  the  logical 
choices  for  representing  the  correction  vector  potential  components  in  the 
far-field  region  (k^r  >>  10).  Working  in  the  spherical  coordinate  system 
(r,6,^),  and  neglecting  all  terms  containing  r ",  r \ •••,  one  can  easily 
show 


7V. ? » -k'll^ r , (5.111 

where  3 is  the  total  vector  potential.  The  total  electr  ic  t ie  Ld  E , i* rom 
Equation  (2.8),  can  therefore  be  shown  to  have  no  r-component,  namelv, 

2 ■>  * 

E(r,3,fi)  =•  (k~  + 77. )3  * k~[(cos  6 cos  ) 1^  - sin  ' 2^)8  - sin  t n $1 

(5.  U) 


As  expected,  the  above  far-field  expression  represents  two  plane  waves 


•Til* 


ante  of  an  unloaded  dipole  antenna  (2L  = 10m)  radiating  in  free  space 


(polarized  in  the  d-  and  the  f-  directions'  propagating  away  from  the  (x,y,z) 
origin  defined  in  Figure  IS.  In  summary,  the  t'ar-electric-f ield  radiation 
pattern  due  to  a current  element  with  no  v-component  can  be  readilv  computed 
by  initially  evaluating  the  total  vector  potential  components  via  the  ROM 
approximation  and  then  using  Equation  (5.1-*'  to  obtain  the  F-field  values 
at  the  desired  observation  points. 

The  far- fie  Id  pattern  for  a given  antenna  structure  is  simply  obtained 
bv  applying  the  superposition  theorem  to  the  individual  radiation  patterns 
of  the  antenna  current  segments  defined  in  the  method  of  moments  approximation 
Radiation  pattern  examples  are  included  for  the  various  antenna  structures 
analyzed  in  the  following  sections. 

5.  -t  Horizontal  Antenna  over  Lossy  Halt-Space 

The  general  developments  of  Sections  5.1  - 5.3  are  applied  to  the 
horizontal  antenna  shown  in  Figure  Jl.  Fortunately,  because  of  the  symmetries 
present  in  this  geometry,  the  [Tlm^l  matrix  in  Equation  IS.*"  takes  the 
following  form  (.Toeplitz  matrix': 


a , 


a , a 


3 


Vi! 


. imp 


O • l M 


n-  l * n- . 


Therefore,  one  needs  to  compute  only  one  row  of  this  matrix  and  use  the  afore 


mentioned  symmetry  to  complete  it. 


The  main  program,  HORtT.  i see  Appendix  IV' 


Medium  I (air) 
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e2r=  <g 
cr2  8 a 


Medium  H (ground) 


Figure  Center- ted  horizontal  dipole  over  a lossy  halt-spaeo. 
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is  developed  to  analyze  the  horizontal  antenna  of  Figure  21.  Using  this 
program.  Figures  22  and  23  are  generated  to  show  the  impedance  variations 
of  a 2L  = 10  meters  center-fed  horizontal  antenna  located  h = 3 meters 
above  various  lossy  grounds.  Radiation  pattern  of  this  antenna  at  15  MHz 
is  also  shown  in  Figures  24  and  25. 

5.5  Vertical  Antenna  over  Lossy  Half-Space 

The  vertical  dipole  shown  in  Figure  26  is  the  next  geometry  considered. 
Unfortunately,  as  was  the  case  for  the  horizontal  antenna,  the  total  [Zlmp] 
matrix  is  not  in  a Toeplitz  form.  However,  main  program  VERT  (in  Appendix  XV) 
is  designed  to  cake  maximum  advantage  of  the  available  symmetry.  As  an 
example.  Figure  27  is  included  to  show  the  radiation  pattern  of  a 2L  = 10 
meters,  h = 8 meters,  center-fed  vertical  dipole  at  resonance  (f  = 15  MHz) 
located  over  various  lossy  grounds. 

5.0  Inverted  Vee-Dipole 

As  a complicated  example,  the  inverted  Vee-dipole  of  Figure  28  is 
considered.  Again,  as  in  the  two  previous  sections,  symmetry  is  used  in  the 
main  program  VEEDIP  (Appendix  IV)  in  constructing  the  [ Z 1IP‘P  ] matrix.  The 
program  is  tested  for  an  inverted  Vee-dipole  structure  having  L = 7.5  meters, 
h = 10  meters,  and  ¥ = 90°;  Figures  29  and  30  demonstrate  the  radiation 
pattern  of  this  structure  at  10  MHz  and  for  various  lossy  grounds. 

In  all  three  of  these  examples,  care  has  been  taken  not  to  violate  the 
conditions  > 10  and  Equation  (4.26)  to  ensure  the  accuracy  of  the  results. 
Also,  since  the  [Zlmp]  matrix  for  these  examples  turns  out  to  be  symmetric, 
a special  inversion  routine  (XINVZ  in  Appendix  IV)  is  employed  to  save  an 
appreciable  amount  of  computer  time. 
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CONCLUSIONS 
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Based  on  the  steepest  descent  path  (SDP)  integration  technique,  an 
efficient  numerical  integration  procedure  is  developed  in  Chapter  3 for 
computing  the  Sommer f eld  infinite  integrals  present  in  the  vector  potential 
expressions  of  a current  element  radiating  over  a lossy  half-space.  Even 
though  this  procedure  is  about  an  order  of  magnitude  faster  than  the  latest 
reported  Sommerfeld  integration  techniques,  the  computation  time  for  a tvpical 
antenna  problem  can  still  become  prohibitive.  The  reflection  coefficient 
method  tRCM)  approximations,  which  are  simply  the  first  term  in  the  asymptotic 
expansion  of  the  Sommerfeld  integrals,  offer  a simple  closed-form  solution 
valid  only  at  the  high  end  of  the  frequency  spectrum  and  which  cannot  be 
employed  in  many  practical  situations.  Also,  the  addition  of  the  second  terra 
in  the  aforementioned  asymptotic  expansion  to  the  RC>!  approximations  is  ruled 
out,  since  the  resulting  vertical  vector  potential  components  diverge  from 
their  respective  exact  integration  values. 

Chapter  4 presents  a novel  approach  in  which  the  transform  domain 
representation  of  the  vector  potentials  is  approximated  such  that  the 
resulting  space-domain  expressions  do  not  require  any  kind  of  infinite 
integration.  This  approach  has  the  merit  of  being  computationally  over  an 
order  of  magnitude  faster  than  the  SDP  technique  of  Chapter  3,  while  being 
accurate  over  a wide  range  of  parameters  of  practical  interest  and,  in 
addition,  offers  a simple  and  numerically  manageable  procedure  for  obtaining 
the  near  E-  and  H- fie  id  components. 

The  general  computer  program,  listed  in  Appendix  IV,  is  developed  hv 
employing  the  approximate  formulas  of  Chapter  -t  and  is  used  to  solve  several 


30 


antenna  geometries.  With  minor  modifications,  this  program  can  be  adapted 
to  analvze  most  three-dimensional  thin-wire  antenna  structures  over  a lossy 


half-space. 
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APPENDIX  I 


EVALUATION  OF  AND  ^ AT  9,  - 0 


In  this  appendix,  the  behavior  of  Equations  (2.47)  - (2.49)  is  studied 
at  9.,  = 0.  In  their  present  forms,  these  integrals  are  not  defined  at 

= 0,  although  it  is  clear  that  their  equivalent  forms  in  (2.42)  - (2.44) 
are  bounded.  Equation  (2.47)  can  be  expressed  here  for  convenience  as 


0 viz  4iTj  J p r-z, 

kA~  - A + /<k“  - A 


(2) 

- Hj“;(p2A)e 


-jz-,  A:-a' 


dA  (I.l) 


where  p9  = r;  sin  (3,,.  The  Hankel  function  in  (I.l)  is  not  bounded  at 

(2) 

0 = 0.  To  circumvent  this  difficulty,  one  replaces  Hq  with  its  expansion 

from  [22] 

To  YP,A  2 r (-l)m+1  [D^Xl2m 

»0  ¥ • V°2X>  - J ?l"  -f  VP2X)  + j l p-T  *<"> 

m=l  (m!)  1 J 


where  y is  Euler's  constant  and  P(m)  represents  the  harmonic  series. 


4>(m)  = 1 + 1/2  + 1/3  + ...  + 1/m  . (1.3' 
Note  that  both  and  the  summation  terms  in  (1.2)  are  even  functions  of 
A,  hence,  their  contributions  to  the  integral  (I.l)  are  zero. 

Substituting  (1.2)  into  (I.l),  one  finally  arrives  at 


.-i  r 

0 viz  4iTj  J 


— » </k~  - A + »K'k“  - A‘ 


- JQ(r,  sin  e2  A) 


• Ln (r .,A)  e 


-jz2/k“-A* 


dA 
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APPENDIX  II 


ASYMPTOTIC  EVALUATION 


In  this  appendix,  a general  formulation  is  developed  for  a higher-order 


asymptotic  evaluation  of  an  integral  with  the  following  format 


u - JL  j p(o  e'jkr  cos(C-0)  d5  , 


(ii. i) 


where  it  is  assumed  that  kr  is  a large  parameter,  — tt / 2 < 0 < it/2,  P(5) 
is  a slowly  varying  function  and  path  f is  shown  in  Fig.  2.  For  large 


values  of  kr,  one  is  usually  interested  in  determining  the  asymptotic 
expression  of  (II. 1);  this  is  done  by  employing  the  method  of  the  steepes:- 
descent  path  integration.  At  the  saddle  point  5 = 9,  one  can  deform  the 
integration  path  F to  the  steepest  descent  path  (SDP)  defined  by 
Re[cos  (5  - 9)]  =1.  Assuming  that  in  this  deformation  no  poles  or  branch 
points  are  encountered,  one  may  express  (II. 1)  as 


P(5)  e‘jkr  Cos(?'6)  d5  . 


(II. 2) 


Since  on  the  SDP  the  relation  Re[cos  (5-9)]  = 1 holds,  one  can 
introduce  the  change  of  variable 


cos  (5  - 02)  = 1 - jt~ 


(II. 3) 


or  equivalently, 


t - n e-J’/4 


(II.  I) 


in  which  t is  a real  variable  taking  the  domain  [-*>,<»].  Substituting  (II. 


into  (II. 2),  one  arrives  at 
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^-jkr-ju/4  r _kcc2 

— j 0(t)  e dt 

2/2  7T  J 


(II. 5) 


Q(t)  = P(£)  sec  — - 


(II. 6a) 


in  which  £ is  replaced  with 


TT  2 /2 

t,  = ± — + jLn(t^  + j + 1 1 >t“  + 2 j ) 


+ 9 , t > 0 


(II. 6b) 


and  Ln  is  interpreted  as  being  its  principal  value.  The  complete  asymptotic 
expansion  procedure  [23]  is  now  used  for  the  asymptotic  evaluation  of  (II.  5). 
In  this  procedure,  one  first  expands  Q(t)  in  a Taylor  series 


Q(t)  . y Q _(P)  tn 

n=0  r(n  + 


(II. 7) 


where  Qv  ' (0)  = Q(t) j and  T is  the  Gamma  function.  Then  (II. 7)  is 

3tn  1 t=0 

substituted  into  (II. 5)  to  finally  result  in 


-jkr-jTT/4  » 0-2n  _n_1/2  (2n) 

l =—r  (kr)  ' QU  ; (0) 

2/2tt  n=0  n' 


(II. 3) 


In  constructing  the  preceding  equation,  the  following  identity  was  used. 


r n -krt2  . 

t e dt  = 


(kr)_(1+n)/2  F[(l  + n) / 2 ] for  n even 


for  n odd 


(II. 9) 


(2n) 

The  task  is  now  to  determine  Q 's  in  terms  of  P.  This  is  achieved  by 


differentiating  (II. 6a)  and  arriving  at 
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r 


Q(0)  = P(9) 

(2) 


V. 


(0)  = 2j 


q(4)(0)  = -4 


+ 1 PJ 

■ dC  J 


5-0 


d4P(C)  5 d2P(C)  , _9_  p 

LdC4  2 d52  16 


»] 


(II. 10) 


C*9 


The  higher-order  terms  can  also  be  determined  in  the  same  fashion.  It  is 
worth  emphasizing  here  that  in  deriving  the  preceding  equations,  the 
following  assumption  has  been  made:  neither  the  poles  nor  the  branch 
points  of  P(£)  are  intercepted  by  the  path  SDP. 

To  present  an  example,  the  higher-order  asymptotic  expansion  of  the 
following  Hankel  function  of  the  second  kind  and  order  v is  derived: 


H(2)tf)  - A f e~JV^w/2  e'jn  COS?  d£  . 
V IT  J 


(II. ID 


Comparing  (II. 11)  with  (II. 1),  one  obtains 

p(0  = 4j  e_jV?+jX;iT/2  , (11.12). 

where  it  is  assumed  that  Q » V.  Substituting  (11.12)  into  (II. 10)  and 
simplifying  the  result,  one  finally  arrives  at 

j Vir/2 


Q(0)  = 4j  e 

-8 


Q(2)(0) 


»24 


jvrr/2 


Q<‘>(0,  . | V2  + £)  e^'2 


(11.13) 


(2) 

The  asymptotic  expansion  of  H is  then  determined  using  (II. 8)  to  be 


V / 1TM 


? + _ ° • ^2(,/  ■ 2V' + &) 


_7/"> 

+ 0(fl  ' ). 

(11.14) 


The  purpose  of  this  appendix  has  been  to  formulate  the  necessary  steps 


for  determining  the  asymptotic  expansion  of  the  vector  potential  components 
expressed  in  (2.47)  - (2.49).  Rearranging  these  three  expressions  into 
the  form  of  Eq.  (II. 1),  one  can  easily  define  a P function  corresponding 
to  each  of  the  vector  potential  components  as 


Pvl2(^  * IvOkIIC 


,(2) 


■ — V (k1P2  sin£)  expfjk^p,,  sint") 
k cos  E,  + /<  - sin'  £ 


(11.15) 


Phlx(J:)  " rhOkl 


sin  E cos  £ ..(2),. 

• — HQ  (k^,,  sin  O exp (jk^p0  sin  O 

cos  t,  + r<  - sin'  £ 


(11.16) 


r 


,,,  ...  . 2 . , cos  t.  - /<  - sin"  £ ..(2)  .. 

Phlz(C)  “ “ j Ih0kl  ' os  sln  ^ cos  ^ ~ Ht  (k1P2  sin  t 


< cos  £ + ►£"  - sin'  £ 


exp  (Jk  p,  sin  £) 


(11.17) 


where  functions  P , , P,  . , and  P,  , correspond  to  the  -IT  , , _ FI,  , , and 
viz  hlx  hlz  0 viz  0 hlx 

..IT,  . expressions,  respectivelv.  Since  the  final  results  take  a 
0 hlz  r 

complicated  form,  the  following  notations  are  introduced  for  the  ease  of 
representation,  namely. 


c =»  cos  0., 


(11.18a) 


s » sin  8., 


q = vK  - sin'  8,, 


(II. 18b) 


(11.18c) 


Performing  the  rather  tedious  differentiation  needed  in  (II. 10)  and 
using  the  asymptotic  expansion  for  the  Hankel  functions  derived  in  (11.14), 
one  finally  arrives  at  the  following  two-term  asymptotic  expression  for 
the  aforementioned  vector  potential  components: 
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“jkir2 


0 viz  vO  <c  + q 4irr„ 


+ XvO  JKC 


A1(3c  - 2/c)  + A.,  + <c 


Kc  + q 


(kc  + q ) " 


2 2 *") 
2S  Aj 


(<c  + q)' 


■jkir2 


AlTkir2 


— + 0(k1ro) 


2c  e 


_jkir2 


0“hlx  ‘hO  c + q 4irr0  + XhO  jc 


-3 


S + A2  * C 


(11.19) 


c + q , , 

(c  + q) 


~jkir2 


7— r + 0(kir2>'3 

4nk^r7 


(11.20) 


-jkir2 


0‘\lz  * rhO'C°S  ^2  SC  Kc  + q 4TTr0 


+ Ih0  cos  <P2  jsc(c  - q) 


<c  + q 


A (5c  - 2/c  + 2s“/q)  + A7  + <c  2s~A“ 

+ — 5 + - 


(KC  + q)' 


(KC  + q)' 


_jkir2 


4iTk^r  “ 


+ O^r,,) 


-3 


wheie  A - A expressions  take  the  following  form: 
1 4 

* Kq  + c 

4 0 3 

A7  = (s  - 2ks~  + K)/q 


A ^ = 2s"/q“  + (3c  - 2/c) /q  - 2 


(11.21) 

(II. 22a) 
(II. 22b) 
(II. 22c) 


A,  = s^c/q^  + s^/q~  + 2(3c  - l/c)/q  - 6 ; (II.22d) 


s,c,  and  q were  defined  in  (11.18).  It  is  interesting  to  point  out  that 
in  deriving  (11.19)  - (11.21),  one  initially  encounters  singular  terms 
at  = 0.  As  expected,  after  some  algebraic  manipulations,  these  singular 
terms  completely  cancel  out  in  all  three  cases,  and  the  two-term  asymptotic 
expressions  for  the  vector  potentials  remain  bounded  for  0 < < 90. 
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APPENDIX  III 

VARIOUS  PARTIAL  DERIVATIVES  OF  g 


The  tree-space  llreen's  function  g,  which  is  defined  as 


2,1/2 


g - gfx^  x.,,  x.?)  - exp  (-Jkr)/4irr  ; r - (x“  + x“  + x“P'“  (III.l) 

is  singular  only  at  the  point  r » 0,  therefore,  its  partial  derivatives 
exist  for  all  r > 0 and  can  be  represented  in  i closed  form.  In  order 
to  express  these  derivatives  in  an  organized  manner  and  suitable  for 


numerical  evaluation,  the  auxiliary  function  R.  Is  introduced  as 


Rj  - (-1)  (I 


-2i 


l - 1,2.3,. 


. ( 2i  - 2)  + jkr  [1  x 3 x 5 x...  x(2i  - 3)]lr‘ 

(III. 2) 


where  R^  has  the  following  convenient  property 


3*.  Ri  “ XJKl+l 


i - 1,2 


1.2.3 


(III. 31 


The  various  partial  derivatives  of  g can  now  be  expressed  in  terms  of 
x.,  R,,  and  lower-order  partials  of  g. 

a)  Partial  derivatives  ot  one  variable 
3 

r — g - x R g 
3x 11 


( 1 1 1 . 4) 


(III. 5) 


~ j g “ ( 3x ( R , + x ? R j ) g + 2(Rj  + x‘R2)  ~ g + x jRj  ~ g 


(.III  . t,  3 


3x: 


where  j * 1.2,3  holds  for  Equations  (III. 4)  - (111-0). 


SO 


b)  Mixed  partial  derivatives  of  two  variables 


* - XJ\R2*  + XkRl  3^7  8 ’ XjXk(Ri  + R2)s 

J k J 


(HI. 7) 


3 ,2 

7i—  8 a Vr2  + x?3)8  + 2xj\R2  dr 8 + xkRl  71  8 (III*8) 

3xj  3xk  J 3xj 


•■n  \ J 


--3-—  g - xk(3xjR3  + xjR,  )g  + 3xk(R2  + *JH3)  3^  8 + Sxjx^,  g 


j k 


33 

+ xkRl  3^  8 


(III. 9) 


— t s ” (Rt  + u:  + xr)R  + x;x"R,]g  + x . (R,  + x-R.)  r~  g 


. - - - - 'j  ' “k'“3  ' “j~ k“4,e  ' k 2 j*'3'  3x, 

3x , 3x,  j j j k 

j k 


* 2.  <8,  * *;«3>  « * w,  * *R»2>  72  * * -*AR2  OT« 

J J 3x.  j k 


+ X,  R 


k 1 , 2 . 
3Hj  3xfc 


(I II. 10) 


where  j i k and  j,k  * 1,2,3  is  assumed  for  Equations  (.111.7)  - (III. 10). 
c)  Mixed  partial  derivatives  of  three  variables 

,3  , 1 


7 - 7 — g = s.x  x,R  j + 2x  x,R,  -r — g + x„R  ~ r — g (III.  11) 

3x , 3x,  3x,  * k 5,  3 k l 2 3x  . i l 3x,  3x, 

' k i ’ 1 >■ 


j 


'i  k 


3" 


3x7  3x,  3x(, 
k e, 


g = x,  x , (R , + x;R, )g  + 3x  x x ,R  v — g + 2x  x ,R , - — 7 g 
k v 3 Jo  J k v 3 3x  k . - 3x- 


+ X X,R,  t -r g + X.R 

J , 2 )Xj  3xk 


l 


^x7  )x. 


(III. 12) 





mm 
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APPENDIX  IV 


COMPLETE  COMPl'TER  LIST  I NO 


In  this  appendix,  a complete  listing  of  a Fortran  program  is  listed 
for  analyzing  wire  antenna  structures  over  a lossy  half-space,  based  on 
the  developments  of  Chapters  a and  5.  Also  included  are  subroutines 
for  evaluating  the  correction  vector  potential  via  the  SDP  integration 
technique  introduced  in  Chapter  3.  To  avoid  confusion,  care  has  been 
taken  to  use  the  same  symbolsand  names  in  the  program  as  were  introduced 
in  the  text,  and  comment  statements  have  been  included  frequently  to 


describe  the  function  of  each  routine. 
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